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Abstract 

The construction of Dirac observables, that is gauge invariant objects, in General Relativity is tech- 
nically more complicated than in other gauge theories such as the standard model due to its more com- 
plicated gauge group which is closely related to the group of spacetime diffeomorphisms. However, the 
explicit and usually cumbersome expression of Dirac observables in terms of gauge non invariant quantities 
is irrelevant if their Poisson algebra is sufficiently simple. Precisely that can be achieved by employing 
the relational formalism and a specific type of matter proposed originally by Brown and Kuchaf , namely 
pressureless dust fields. Moreover one is able to derive a compact expression for a physical Hamiltonian 
that drives their physical time evolution. The resulting gauge invariant Hamiltonian system is obtained 
by Higgs - ing the dust scalar fields and has an infinite number of conserved charges which force the 
Goldstone bosons to decouple from the evolution. In previous publications we have shown that explicitly 
for cosmological perturbations. In this article we analyse the spherically symmetric sector of the theory 
and it turns out that the solutions are in one-to-one correspondence with the class of Lemaitre-Tolman- 
Bondi metrics. Therefore the theory is capable of properly describing the whole class of gravitational 
experiments that rely on the assumption of spherical symmetry. 
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1 Introduction 



Until today General Relativity (GR) is the most successful theory in describing observations involving 
gravitational interactions on macroscopic scales. It has been tested in various experimental settings and so 
far many predictions of GR have been confirmed to great accuracy. But next to the complicated highly non 
- linear structure of Einstein's equations there is one aspect in GR which has to be treated with special 
care as compared to other field theories such as Maxwell's theory: the issue of observables. Observables are 
those quantities that respect the gauge symmetry of the theory, for Electrodynamics this would simply be 
all U{1) - invariant quantities. In contrast for GR the gauge group is closely related to the group of space 
time diffeomorphisms, denoted by Diff(M), reflecting the independence of physics on space time coordinates. 
We are dealing with a background independent theory where "space" and "time" do not have any a priori 
physical meaning. 

In the canonical picture of GR [1], where one performs a (3+1) - split of the four dimensional space time 
(M, g) , gauge invariance carries over to the condition that phase space functions have to commute with 
the Hamiltonian and spatial diffeomorphism constraints of The reason why it is so complicated to 

extract the gauge-invariant content of GR is that the Poisson algebra formed by the constraints of GR, 
the so called Dirac algebra, is an extremely difficult one. It is not only non-Abelian but is not even an 
honest Lie-algebra: instead of structure constants it involves structure functions, i.e., phase space dependent 
quantities. However, even if observables were known there would still be a conceptual task to solve which 
is often referred to as the problem of time in GR [H |5] and concerns the evolution of observables. Speaking 
about dynamics of GR we usually refer to Einstein's equations. In the canonical picture these can be obtained 
from a canonical Hamiltonian and the Hamiltonian and diffeomorphism constraint equations. However, as a 
consequence of diffeomorphism invariance the canonical Hamiltonian is a linear combination of constraints 
onljl. Hence, the Hamiltonian equations generated by this canonical Hamiltonian will be identical to 
zero in the case of observables. Consequently, Einstein's equations do not describe what we usually would 
call physical evolution but rather describe how the metric changes under gauge transformations. What is 
needed for observables is a gauge invariant version of Einstein's equation, generated by a so called physical 
Hamiltonian and describing non - trivial evolution of observables. 

A framework that allows to construct observables and analyse their evolution for constrained systems such 
as GR is the so called Relational Formalism. The basic idea is to take the background independent nature 
of GR seriously and define observables not with respect to unphysical space time points in M but to use 
relations between dynamical fields instead. These ideas date back to the seminal work of Bergmann and 
Komar [6l[71[8] from the 1960s. Its conceptual foundations were very much improved in the 1990s, see [9l[T0] 
and references therein. The corresponding mathematical framework was also developed in the 1990's, see 
e.g. [11] and references therein and rediscovered more recently in [121113). Once appropriate dynamical fields 
are chosen as clock^ in order to give space and time a physical meaning one can at least formally write down 
an expression for observables associated to any phase space function /. In [HKTS] a perturbative scheme to 
compute these observables was developed and its application to perturbation theory around Minkowski space 
and cosmological perturbation theory for different choices of clocks was discussed. However, the explicit 
form is a power series in these clock variables with coefficients involving multiple Poisson brackets of the 
constraints and / generally leading to a rather complicated dependence on the physical time parameter that 
in most cases cannot be written down in a closed analytic form. 

However, as pointed out for instance in [16], the precise expression of the Dirac observables in terms of 
non gauge inavriant objects, which generically involves a hopelessly complicated, infinite series of multiple 
Poisson brackets, is in fact irrelevant from the point of view of the physical (or reduced) phase space that 
we are actually interested in. What is relevant is the Poisson algebra of those Dirac observables and the 

^See m or [3] for a detailed description how to deal with gauge symmetries in physical theories. 
^In the absence of boundaries such as in asymptotically flat situations. 

^Whenever we talk of clocks we mean devices (fields) to measure temporal and spatial distances. 
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physical Hamiltonian which drives their physical time evolution. General expressions for an arbitary first 
class system and for general choices of clocks were derived in |16j . From a practical point of view and 
especially with regard to quantisation it is of course desirable to have a sufficiently simple gauge invariant 
Poisson algebra Slphys and a compact expression for the resulting physical Hamiltonian Hphys- The structure 
of both 2tphys, Hphys depends sensitively on the choice of clocks and for generic choices both are beyond 
mathematical control. One therefore should analyse first which type of matter simplifies their structure. 

First steps in this direction have been performed in [T^ where a phantom scalar field, i.e. a scalar field 
with negative kinetic energy term, has been used as a clock and the corresponding physical Hamiltonian has 
been derived. The choice of a scalar field with vanishing potential as a clock field was motivated entirely by 
mathematical considerations based on the Brown - Kuchar^ mechanism jl8j which allows to deparametrise 
the Hamiltonian constraint. From a physical point of view, introducing some ad hoc matter component 
that adopts the role of clocks for GR might look artificial and even dangerous, especially if that matter is 
in conflict with the ususal energy conditions that stabilise the system. However, the way that the relational 
formalism works is in fact very similar to the Higgs mechanism: In the gauge invariant formulation, the 
scalar field completely disappears. What remains is the corresponding physical Goldstone boson and what 
one has to worry about is that this additional degree of freedom is not in conflict with observation, in 
particular, that the energy conditions hold in the gauge invariant description. 

A more physically motivated choice of clocks which have the additional advantage of implementing the 
Brown - Kuchaf mechanism is defined by the original Lagrangian for pressure free dust due to Brown 
and Kuchai^ [18]. This choice of matter is a sense physically distinguished because it can be considered 
in a precise sense as a congruence of mutually non interacting, freely falling observers which only interact 
gravitationally and which define the dynamical reference frame for GR. From the mathematical point of 
view, the essential point in the construction is that the constraints of the coupled system, including gravity, 
matter and dust, can be written in a deparametrised forrrfl. As was shown explicitly in |19j . this achieves 
the goal of drastically simplifying the Poisson algebra of physical observables and the physical Hamiltonian 
and furthermore ensures that the associated physical Hamiltonian is time independent. In contrast to |18] . 
the requirement that the physical Hamiltonian should be positive definite, or equivalently physical time is 
required to run forward rather than backwards, demands the dust needs to be phantom dust. This is exactly 
the same reason why the phantom occurred in the previous work [17] and ensures that the energy momentum 
tensor of observable matter obeys the usual energy conditions. It was also shown that the resulting theory 
is in good agreement with current cosmological observations in [20] as far as FRW space times as well as 
perturbations around FRW are considered. The reason for why that happens is the afore mentioned Higgs 
mechanism by which the phantom dust completely disappears from the physical particle spectrum together 
with another nice feature of this particular choice of clocks which is absent for generic choices: The physical 
Hamiltonian system posesses an infinite number of conserved charges which can be considered as physical 
energy and momentum densities of the dust respectively and which enforce the corresponding gravitational 
Goldstone boson modes to decouple in a mathematically precise sense from the physical time evolution. In 
particular, it is consistent to tune those charges to be arbitarily small so that their corrections to Einstein's 
equations for the non Goldstone modes are also arbitrarily small. In this way the dust comes as close as 
possible to the mathematical idealisation of a test observer while taking its gravitational interaction into 
account. Of course, one could also use gravitational clocks, that is certain components of the metric tensor 
and thus avoid the Goldstone modes altogether. However, while that works well for the linearised theory, 
taking into account the full non linearities of Einstein's theory leads to equations of motion for the physical 
graviational modes that have no resemblance with the classical Einstein equations whatsoever. Similar to the 

*A system of constraints Ci is said to be deparametrisable if one can find a local coordinate chart on phase space with 
two mutually commuting sets of canonical pairs denoted by {q°',Pa) and (T^ ,Pi) such that in this chart the constraints can be 
written in the locally equivalent form Ci = Pi + hi where hi depends only on {q°',pa)- Equivalence has to be understood in 
the sense that the constraint hypersurface, spanned by the deparametrised constraints is exactly the same as the one spanned 
by the original set of constraints. 
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vast number of proposals from cosmologists, particle physicists and string theorists for possible extensions of 
the field content of the standard model such as inflatons, axions, dark matter, supersymmetric extensions, 
dilatons, Kaluza Klein modes etc. we consider the existence of a fundamental dust field in addition to 
standard matter as an at least theoretical possibility which of course has to be tested in experiments. The 
fact that the dust only interacts gravitationally makes it in principle a perfect dark matter candidat^. 

The topic addressed in this article is the spherically symmetric sector of the theory presented in |19j . 
We will show in detail that a gauge invariant version of Einstein's equations specialised to the spherically 
symmetric case can be mapped to a family of Lemaitre - Tolman - Bondi (LTB) solutions. The only 
modification is that the equations contain a phantom dust energy density instead of the usual dust energy 
density and we will discuss explicitly which consequences this has. 

As pointed out in [16], the existence of gauge invariant version of Einstein's equations and thus a re- 
duced phase space formulation of GR involving an algebra of observables and a physical Hamiltonian is 
of advantage when a quantisation of gravity is concerned because one sidesteps the difficulties involved in 
the anomaly free quantisation of constraints and the construction of the physical Hilbert spaces. Recently, 
this strategy has been used in order to present a reduced phase space quantisation [21] for Loop Quantum 
Gravity [22^123]. This framework developed for the full theory can of course be applied also to the spherically 
symmetric sector |24j . 

The paper is structured as follows: 

In section [2] we review the results of [19] for the full theory. In section [3] we specialise the observables 
and the physical Hamiltonian to spherically symmetric space times and in section [4] derive the correspond- 
ing equations of motion. These equations are a gauge invariant version of Einstein's equations in the case 
of spherical symmetry. Section [5] discusses the solution of these equations which belong to a family of LTB 
- solutions. In section [6] and [7] we discuss further properties of these solutions such as their semistatical 
behaviour as well as the occurrence of singularities. For the latter we analyse in particular the effect of 
the phantom dust energy momentum that occurs as a source in these equations. Finally in section [8] we 
conclude, discuss the implications of the work done in this paper and give an outlook. Several appendices 
on the boundary conditions imposed on gauge invariant variables, various spherically symmetric coordinate 
systems employed in the literature and a comparison with the covariant derivation of the LTB solutions 
complete the paper. 

2 Brown — Kuchaf Dust Reduction of General Relativity 

In this section we summarise the analysis performed in [19j which builds on the seminal work |18j . 
The Brown - Kuchaf dust Lagrangian is given by 



^One could call it a NIMP (Non Interacting Massless Particle) as compared to a WIMP (Weakly Interacting Massive 
Particle) which is considered as the most favourite (cold) dark matter candidate. Massless here refers to the fact that the dust 
Lagrangian does not contain a usual mass (or potential) term. However, the foliation defined by the dust is always spacelike 
and the corresponding foliation vector field has unit timelike norm with respect to the physical metric (but it is not necessarily 
normal to the foliation). That is, the dust moves along unit timelike geodesies which also follows from the fact that the dust 
Lagrangian can be interpreted as the sum (or integral) of relativistic point particle Lagrangians (one for each flow line of the 
congruence) with variable mass depending on the flow line. This mass distribution is just the energy density of the dust. See 
[18j for details. Hence the dust moves at non relativistic speeds as required by realistic dark matter models. Of course the 
dust energy density has the wrong sign to explain the anomalous galactic rotation curves but we should recall that the dust 
is anyway like a Higgs boson which disappears from the observable particle spectrum. There is no obstacle in adding one's 
favourite observable dark matter model. 
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which is coupled to the standard Einstein Hilbert action Seh 



Seh = - [ d^XvlditM^'^^- (2.2) 



K 
M 

Standard matter Lagrangians can be added in the usual way, but they do not couple to the dust. Here M 
is a four dimensional manifold which can topologically be identified with M x A' for some three dimensional 
manifold X of arbitrary topology and X ^ M are local coordinates on M. Next, g^y{X) with = 
0, 1, 2, 3 denotes a (Pseudo-) Riemannian metric on M. U £ T*M is a one form defined as the differential 
U = -dT + WjdS^, j = 1,2,3, for some scalar fields T,Wj,S^ G C°°(M). So finally the action is a 
functional of (j^j/ and eight scalar fields p, T, Wi, 5*. ^^^R is the Ricci scalar corresponding to the metric g. 

We need the Hamiltonian formulation of that system which can be derived by performing a usual ADM 
- (3+1) - split with respect to a foliation oi M = M. x X. We denote the momentum conjugate to the 
ADM 3 - metric Qab by p""^. The momenta conjugate to the configuration variables p, Wj,T, are denoted 
by Z, , P, Pj respectively. Latin indices range from 1 to 3. A detailed Dirac analysis of the occurring 
constraints shows that the coupled system is second class. Hence, in order to proceed, one passes on to the 
corresponding Dirac bracket and solves the second class constraints explicitly. As a consequence it turns out 
that the momenta conjugate to p and Wj vanish and that these configuration variables can be expressed in 
terms of the remaining phase space variables. Explicitly, we have 



Z := 0, Z^ := 0, Wj := p' := -== q'^'UaUb + 1 . (2.3) 

P V det q L J 

In general, analysing second class constraint systems is a very hard task due to the complicated structure of 
the corresponding Dirac bracket. However, it turns out that for the system Seh + Sdust, when restricting 
attention to the geometry variables {qab, p''^) as well as the remaining dust variables {S^,Pj) and (T, P), 
then the Dirac bracket in this sector reduces to the standard Poisson bracket again, which makes a further 
analysis of this system tractable. 

We end up with a first class system possessing the following constraints 

gtot ^ c + c^"*^* (2.4) 
cr = ca + ct^' , (2.5) 

where the geometry and dust contributions can be written as 



c 

^dust 
„dust 



--qacDbp'"' (2.6) 

-^i= [qacqbd - \qabqcdYV " Vdet(g)i?(<?) (2.7) 
K A/det q 2, 

PyJl + qabUaUb (2.8) 

PTa + P,S^,a ■ (2.9) 



Here Da is the covariant differential compatible with qab and R{q) denotes the Ricci scalar of qab- Here P 
takes only non positive values and thus the energy density of the dust is negative or zero which is why we 
call it phantom dust in contrast to |18j . The reasons for that can be summarised as follows: If P would take 
positive values then on the constraint surface we would have c < 0. As we will see, the derived physical 
Hamiltonian is approximated by |c| = — c which in the limit of flat space would be the negative of the 
standard model Hamiltonian. One could cure this by letting time run backwards and defining the physical 
Hamiltonian by — |c| = c but then energy would be unbounded from below. Furthermore, the dynamical 
foliation generated by the dust would be past oriented if we chose the other sign, see [19] for a detailed 
discussion. 
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Concerning the interpretation of (j2.ip we refer the reader to [18], where the authors describe that the dust 
action can actually be derived as a field theoretic generalisation of the concept of free massive relativistic 
particles moving on geodesies of the gravitational field created by the entire collection of particles. To get 
an idea why this interpretation is possible, one can check that the integral curves of = g^^Uy describe 
geodesies of 5^1,, is constant along each of these curves and T describes proper time. So = const. 
labels a geodesic and T = const, is an affine parameter along the geodesic. This is exactly the reason why 
it so convenient to use the fields T, as a physical reference frame. 

The canonical Hamiltonian that generates Hamiltonian equations for {qah^p""^)-, {P-, T) which are, together 
with the constraints in (12. 4p and (12. 5p . equivalent to the 10 Einstein equations obtained in the Lagrangian 
framework is given by 

ll^an-= j d^x{n{x)c{x)+n^{x)ca{x)) . (2.10) 
X 

where n is the so called lapse function and the so called shift vector which play the role of Lagrange 
multipliers. 

Now we come to the crucial observation made by Brown and Kuchaf in [18] that the system Seh + Sd is 
indeed (partially) deparametrisable. To see this, note first, that the total Hamiltonian constraint (|2.4p can 
be written in equivalent form 

/ „a6^dust «dust j nO,b^ n, 

c*°* = c + pJl+ '^ ""^Z" ^c + P^l + ^-^ . (2.11) 

On the constraint hypersurface the constraints c*°* and cjf* can be solved for the momenta P and Pj 
respectively and thus be written down in completely equivalenl[£| form as 

-tot ^ h=y'c^- q'^^CaCi, (2.12) 

cf' = Pj+hj = S'}[Ca - hTa] , (2.13) 

where we assumed that S'^a is non degenerate and defined its inverse Sj. As shown in [19] this condition is 
gauge invariant under the gauge transformations generated by (j2.12p and (|2.13p . 

At least the total Hamiltonian constraint is in deparametrised form now, because h depends only on the 
gravitational variables qah^p""^- Moreover, since the constraints are linear in the dust momenta, (j2.12p and 
(|2.13p form a strongly Abelian first class constraint algebra. As a consequence also {h[x), h{y)} = because 
h commutes with P. However, only the total Hamiltonian constraint c*°* is of deparametrised form but not 
the total diffeomorphism constraint, so Poisson brackets between either h{x) and hj{y) or hj{x) and hj{y) 
will not vanish in general. So we achieved a partially deparametrised form for the coupled system of gravity 
and dust. Note that this can also be obtained when additional, for instance standard model, matter would 
be coupled to gravity and dust. The only difference in this case will be that the constraints c and Ca above 
will then not only include gravitational contributions but consist of a sum of gravity plus the additional 
matter contributions. An example where additional to gravity and dust a K.G. - scalar field was considered 
can be found in |19^ I20j. 

We are now in the position to define distinguished coordinates on the reduced phase space defined by 
the first class system defined by (|2.12p and (j2.13p . It is clear that T, are pure gauge and that P, Pj can 
be solved for the gravitational (and standard matter) field variables. Therefore it is natural to introduce a 
four parameter family of gauge fixing conditions defined by T{x,t) = r, S^{x,t) = where r G M and a 

^By equivalence we mean that the constraint hypersurface generated by both constraints are the same. 
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takes values in the dust space S = S{X) which by assumption is diffeomorphic to X. The dust space labels 
the geodesies and r is an affine parameter along the geodesies. 

The relational framework can now be applied and the results can be described as follows (see |19| for all 
the details and the derivations): 

Let Sj{x) be the inverse of S-^aix) and J{x) := det{dS/dx). Consider 



qij{a) = J d'x\Jix)\5ia^ -S^{x))SfS'jqab 

X 

fi{a) = (fx\J{x)\5{<Ti - S^{x)) ' y (x) 

X 

(2.14) 

It maybe checked explicitly that 

{f^ia),m{^')} = K5l,6i^5{a,a') . (2.15) 

The interpretation of (|2.14p is obvious: These are the coordinate transformations of {qahiP""^) respectively 
into the dynamical coordinate system defined by Xg- = S~^{a). One can arrive at these expressions inde- 
pendently by symplectic reduction which explains why {qij,p^^) continue to be a conjugate pair. 
Next let 



°° 1 

Qjk{T,cr) := —{hT,qjk}{n) 

n=0 
CO ^ 

P^\T,a) := Y,-{h^,p>^]^^^ (2.16) 



n! 

n=0 



where we introduced 



j d^a{T -f{o))h{a) with h = h{qjk,p'^) (2.17) 



It maybe checked explicitly that (|2.16|) has vanishing Poisson brackets with all constraints. 

These expressions can no longer be described in a compact form, they are hopelessly complicated to 
evaluate as functions of q,p. However, as stressed in the introduction, all we need is their Poisson algebra 
and their time evolution. To that end, set Qij{cr) := Qij{T = 0,0"), P'^^{a) := P^^{t = 0,cj). Then it maybe 
checked explicitly that 

{P'^{a),Qki{cT')} = K6l,Si^5{a,a') . (2.18) 
still form a canonical pair, all other Poisson brackets vanishing. Furthermore, let us define 



H{a) = ^C^-Q'^CiC,{a) , (2.19) 

where 

C := c{qjk = Qik^f"" = P'^), Q := Cr{qjk = Qjk,!^'' = P'""), (2.20) 

and c, Ci are just c, Ca expressed in the dust coordinate system. Then it may be checked that the physical 
Hamiltonian generating time evolution for all observables fiQij,P^-^) is given by 

Hphys := J d?aH{a) . (2.21) 
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When looking at its variation 



mphys = jd^a[(^^yC-(^^^yC, + ^CiCjQ"'Q^'6Qij] (2.22) 
s 

=: J d^a(^N6C + N'6Ci + ^HN'N^Qij^ , (2.23) 
s 

one sees that the physical equations of motion generated by Hp^ys are almost equivalent to the ones generated 
by the canonical Hamiltonian Hcan with the identification qab{x) — > Qij{(7),p"'^{x) P^^{a) modulo the 
following important differences: First, lapse and shift A^* are not phase space independent functions as 
for Hcan where they only encode the arbitrariness of the foliation. Rather they are observable phase space 
functions composed out of the elementary fields Qij,P^^ as 

N:=^, N^ = -9!^ , (2.24) 

Second, there is one additional contribution proportional to the Hamiltonian density H{a). But H{a) is a 
conserved quantity in the theory and can be freely chosen on the initial value hypersurface, so we may tune 
this term to alter the equations of motion as little as we like. 

The physical Hamiltonian has an infinite number of conserved charges, namely energy and momentum 
density H(a), Cj{a). The latter ones generate active diffeomorphisms of the dust space 5, they are to 
be considered as symmetries of the system rather than gauge transformations generated by the passive 
diffeomorphisms of X . Likewise, the former are related in an intricate way to time reparametrisation 
invariance in General Relativity. 



3 Spherical symmetry 

Now we want to specialise the general theory to spherically symmetric spacetimes. We will work directly at 
the gauge invariant level and assume that what one usually measures in physical experiments are not the 
gauge variant three metrics qab and their canonically conjugate momenta p"'^ in some unphysical coordinate 
system but rather the physical metrics Qij and their physical canonically conjugate momenta P*-' measured 
with respect to the physical reference frame S given by the dust fields. There is nothing to debate about 
the fact that whenever we perform experiments we measure physical gauge invariant quantities and not 
the kinematical quantities qah^p""^- Hence we will require spherical symmetry with respect to the physical 
coordinate system a. 

Thus, spherically symmetric spacetimes M = R x 5 will be characterised by a triplet of Killing vector 
fields {'^i,'^2)C3} S whose commutator algebra is isomorphic to the Lie algebra so(3). As usual, for such 
spacetimes it is always possible to find a coordinate chart in which the physical four metric G^^ takes the 
special form 



G 



-iV2 + Q.jN^N^ 



Qij 



(3.1) 



where i,j = 1, 2, 3 and 



Qij ('^r) 



diag 

sin (Tgdiag 



PKiOr) PdcTr) Pr(cT,) 



2A(cj^) ' 4R{ar) ' 4R{ar) 



sm 



(3.2) 
(3.3) 
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where (cTj., erg, cr^) are spherical coordinates of a physical coordinate system on the spatial slices and (A, Pa) 
and (R, Pr) respectively are conjugate pairs. Hence, the only non vanishing Poisson brackets are given by 



(3.4) 



In contrast to the gauge variant formalism lapse N and shift are not arbitrary phase space independent 
functions but are given by (j2.24p . 

Using these fields the geometry parts of the Hamiltonian and spatial diffeomorphism constraints reduce to: 



C 



Ccjr 



— sin agAR 

K 

1 

— sin ag 

K 



Pk PkPr 2 /R'Y 4 P" AR'A' 2 
APa' + R'Pr 



(3.5) 

(3.6) 
(3.7) 



Hence, the physical Hamiltonian Hphys, which generates the evolution of observables, specialised to the case 
of spherically symmetric spacetimes reads as 



H 



phys 



(3., 



4 Equations of motion 

Now we want to discuss the physical equations of motion for the spherically symmetric case: Physical time 
evolution is generated by the physical Hamiltonian (I2.2ip . so as a first step we need to compute its variation. 
This was already demonstrated for the general case at the end of section [2] and now we want to specialise this 
result to the spherically symmetric case: For this purpose we must use the physical Hamiltonian Hphys that 
has been specialised to the spherically symmetric case in (|3.8p and compute its variation while considering 
A, R, P\ and Pr as the basic variables. This results in the following first order equations of motion: 

A = {Hphys, A} R = {Hphys, R} 

Pa = {Hphys, Pa} P^, = {Hphys, Pi?} • (4.1) 

The explicit derivation of physical equations of motion from the variation (5Hphys of the physical Hamiltonian 
has been done for full GR in [19] and it was shown that suitable boundary conditions can be chosen so that 
boundary terms can be neglected in that calculation. This carries over to the spherically symmetric case. 
A discussion of boundary conditions in the latter case is given in appendix |Al Denoting derivatives with 
respect to r and ar with dot and slash respectively the first order equations of motion explicitly read as 

A = - + (AT^'-VA + iV'^'-A' (4.2) 

4R^ AR ^ ' ^ ^ 

R = -^^^N'^rR' (43) 

• _ NC{ar) NPaPr AN'RR' ANR"R ANRR'A' 
^ ~ i^^A 4AP A^^^V2 X3 

+Af'"'-PA' - {N^-fAH (4.4) 
. _ 2NC{ar) NAPa^ SNPaPr ANR"^ ANR" ANR'A' 
^ ~ 4^^;^^^P3 lR^^n^^~ A^ 

ANA AN"R AN'R' AN'RA' , „ „ , , 
X — + ^;^ + iN'^')PR + N'''PR ■ (4.5) 
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We obtain the same equations when we speciahse the first order Hamiltonian equations for full GR given 
in [19] to spherical symmetry and additionally set the K.G. - scalar field contributions to zero. Formally 
these equations of motion coincide with the ones generated by Hcan for the gauge variant 3-metrics and 
momenta up to the term proportional to the Hamiltonian density H in (j4.4p . but there is one important 
difference: N and N'^'^ are no free functions anymore. Rather they are phase space dependent functions 
given by (j2.24p and choosing different A'^, A^* is equivalent to working in a different physically distinguishable 
coordinate system. This again is due to the fact that Diff(5) is the group of active dijjeomorphisms on S 
as opposed to the passive diffeomorphism on X. 

The Hamiltonian density H{a) as well as C^^ are constants of motion jT9]. The latter corresponds to the 
momentum density of the dust denoted by e'^^^a), that is C'^^{a) = —€'^''{a). In the following we will restrict 
our discussion to the case of vanishing dust momentum density. Then the shift vector N"^ being proportional 
to C"^'' is vanishing. Furthermore the Hamiltonian density is given hy H = y^C^ — Q^^CiCj = C and we 
automatically get unit lapse = 1. This means that equations (j4.2p - (j4.5p simplify significantly, all terms 
proportional to N^'' , N' and N" vanish. Using that H{a) is a constant of motion, meaning that it does not 
evolve in physical time r, we can write 

H{a) = smage{ar) for e(fJr.) > . (4.6) 



From H{a) = —P{a) = /5dust\/det((5) we obtain 



smo-ge(a^) e{ar) , . 

/?dust(T,o-rj = -===— = -— —J- . (4.7) 

VdetQ A{T,CFr)R^(T,ar) 



Hence, we can write 



C(o"r) := j dagda^C{a) = — J dagda^ sin agA{ar)R'^ {(Jr) poiT, CTr) = —4:TTpo{o'r)J^{o'r)R'^{o'r) ■ (4.8) 

For further analysis we want to transform the first order system ()4.2p - ()4.5p into a second order system. To 
this end one can solve (j4.2p and (j4.3p for the momenta and then use these expressions and their r-derivatives 
in ([O]) and (fi3]) . Setting N = 1 and iV'^'' = and using (fi^ this results in 

A R_i^_i^_(RY i(R^y 2r;^_±_ ^ _p_D 

R AR A3i? \rJ ^aArJ ^A2i? i?2 2 ' ^ ' 

2R 2R" 2AR 2A'R' ^ _pD . . 

R ^ A^R AR A^R 2 ' ^ ' 



Furthermore we can use that H{a) = — Pdust\/det((5) yielding 



_ /i?\2 1 2i?A 2i?'A' 2R" 1 _ pD 

~\r) ^1?\r) ~TR~l^^J?R~m ~ ~T ' ^ ' 

and the vanishing of the momentum density 

R' - ^R' = . (4.12) 
A 

(|4.9p - (j4.1ip are actually not entirely independent of each other, in fact it is sufficient to keep (|4.1Up and 
()4.1ip : One can check that solutions to these two equations also satisfy 
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5 Solving the equations of motion 



In this section we will discuss the solution of the equations of motion. We will mainly follow [25] for the 
derivation of this solution using a slightly different notation here. In order to solve (j4.9|) - (|4.12|) it is first 
important to notice that the (physical r-) time dependencies of A(r, 0"^) and -R(r, 0"^) are not independent 
from each other. We can solve (|4.12|) and get: 

A(r,cT.) = i?'(r,cT.)-=i==, (5.1) 

where E{ar) is a so far arbitrary, time-independent functioE0 of a^. Exploiting this, we can rewrite (j4.10p 
and as 

i? R R E py) 

R ~ Rm ^ 2R'R ~ r ■ ^^'^^ 
First, by multiplying (fOj) by R'^R we get 



2RRR + R^-ER = R{2RR) + R{R^ - E) = ^ \r{R^ - E)] =0 

^^7 = , (5.4) 

where we defined U = R{R^ — E). U does not depend on r so we can write C/(r, ar) = F{cjr) for an so far 
arbitrary cir-dependent function F{(Tr)- By differentiating U with respect to (Jr we get 

U' = F' = R'{R^ -E) + 2RRR' - E'R . (5.5) 

Dividing this by Rj^R' we see that F' is equal to the lefthand side of (15. 2p minus 2 times the lefthand side 
of (15. 3p . Thus, we can rewrite F' as 

F' = ]^R>R^p^ = ]^^/i^KR^pj^ = -]^^/iTEe{ar) , (5.6) 



where we used (j5.ip in the first equality and (14. 7p in the second one. Integrating the equation above once, 
we end up with 

F{ar) = -\j d\e{\)^l + E{\)+a=: -iM(a,) + ^M(O) , (5.7) 



where M(0) is an so far arbitrary (constant) mass. M{ar) can be interpreted as the effective gravitating 
mass (multiplied by 1/4 since the relation between the Schwarzschild radius and the central mass involves 
a factor of 2) of the dust inside a sphere with radial label cr,.. So finally the equation U = F yields 



R 

We arrive at the following general form of the metric 



R = ±\l^ + E . (5.8) 



1 + E[ar) 



/ F(ar) 
R{T,ar) = ±jE{ar)+ ^ ' 



R{t, ar 

I 

R'{T,ar)R'^iT,ar) 



2F'{ar) 

pD{T,ar) = p,/ ^ ^p2/ ^ ^ > (5-9) 



'it will become clear later on why we have chosen this rather complicated form for the integration constant. 
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with E{ar) > — 1 and F{ar) must be chosen such that the expression under the square root in the second 
hne is non-negative for all i?(r, (Jr), see section [7] for further elaborations on this topic. 
The solution for /2(r, 0"^) for general E{ar) can only be given in parametric form, distinguishing the cases 
E>0,E = 0,E<0 and using the conventions of [26] one finds: 

• E{ar) > : 



R{t, ar 
sinh(?/) — r] 



■^cosh(r7) — 1 



2E(ar 



2[E{ar)]H(3iar) 
F{ar) 



(5.10) 



Eiur) < : 



R{t, ar) 
rj — sin(ry) 



F{ar) 



1 — cos(?7) 



2{-E{ar)) 

2[-E{ar)\^{f3{a 
F{ar) 



(5.11) 



E{ar) = : 



R{T,ar) 



F{(3{ar) - t) 



2/3 



(5.12) 



Here P{ar), E[aj.) and F{ar) are so far arbitrary functions of cr^ which allow a coordinate choice by a partic- 
ular form of P{ar) and two further physical quantities. As mentioned above F in (jS.Th can be understood as 
the effective gravitating mass within radius ar and E{ar) determines the time evolution of R as well as the 
local geometry. For a more detailed discussion see for instance |27j . The time r at which i?(r, 0"^) is equal 
to zero is j3{ar) and one calls r > j3{ar) the big bang time whereas r < I3{ar) is referred to as the recollapse 
time. Often one refers to these three cases as elliptic, parabolic and hyperbolic for E{ar) < 0,E{ar) = 
and E{ar) > respectively. The solution in (j5.9p can easily be identified as the whole class of Lemaitre- 
Tolman-Bondi (LTB) metrics (see appendix lB.3p . so we managed to map the problem of solving the physical 
r-evolution of gauge-invariant three-metrics Qij to the problem of analysing the LTB class of solutions in 
standard GR. However, it is important to note that opposed to the standard LTB models, the mass term 
M{ar) enters with the opposite sign in F{ar) = |(M(0) — M{ar))- That means the effective gravitating 
mass will decrease when going further away from the central mass. This happens because the dust clocks, 
which have been chosen as a physical reference system, have negative energy. 



We can fix the coordinate choice included in (3{ar) and the mass M(0) in ()5.7p by requiring that in the 
vacuum case the general solution above should reduce to the Lemaitre solution [28] given by 



-dr^ + 



-RLcm(T, ar)) 



dar"^ + Rl^^^{T,ar)d^^ with i?Lem(T, cr^ 



2 V -Rs(fr 



2/3 



(5.13) 



where Rg = 2MG/(P' denotes the Schwarzschild radius which simplifies in units h = c = G = l to Rs = 2M. 
The Lemaitre solution is diffeomorphic to the Schwarzschild solution and can be obtained from the latter 
through a coordinate transformation into a comoving coordinate system. So this metric describes the local 
coordinate system of a free falling observer in a spherically symmetric gravitational field originating from a 



12 



central mass M, more details can be found in appendix IB. 1[ Matching the Lemaitre with the LTB solution 
we obviously need to consider the case E{ar) = 0. This is often referred to as the marginally bound case. 
Furthermore, in the vacuum case, e((Tj.) is zero and as a consequence F' vanishes. Hence, in this case F is just 



and identifying M(0) = 4M because = 2M = 1/2M(0). 
5.1 Newtonian limit 

Let us have a look at the Newtonian limit of this spacetime: In [25] it was argued that i?(r, ar) is the quantity 
that should actually be identified with "Newtonian distance", for instance in the context of luminosity 
distance. Having this in mind let us look at equation (j5.8p again. By differentiating this one with respect 
to r we obtain 



which coincides formally with Newton's equation of motion for a point particle under the influence of a 
central mass. The crucial difference is that the effective gravitating mass F/2 is less than the central mass 
M = Rs/2. So the gravitational field is weaker than one would expect without having the dust around. 
There is one remark necessary concerning the falloff behaviour of R{t, ar)'- One might wonder whether the 
boundary conditions discussed in appendix |A] are not violated by the explicit form of R as in (j5.12p . The 
solution to this puzzle lies in the fact that in the derivation of the falloff conditions ()A.3P we assumed a 
coordinate system which approaches a flat spherical Minkowskian one in the asymptotic limit. In contrast, 
the coordinate system we are using here is the coordinate system of a freely falling observer so one cannot 
expect (|5.12p to hold automatically. Nevertheless, we assumed that e((T,.) — > for — > oo, so in the limit 
of vanishing dust density and in the marginally bound case (that is E=0) the LTB coordinates merge into 
Lemaitre coordinates which themselves can be transformed to Schwarzschild coordinates having the correct 
Newtonian limit (see appendix IB. ip . 

We have seen above that the exact Schwarzschild solution is obtained for F =const. and E = How- 
ever, it is possible to obtain the Schwarzschild solution also for certain E ^ and F =const. The idea is to 
transform the LTB coordinates first into generalised Painleve - Gullstrand coordinates (GPG), introduced 
in [29] and then to Schwarzschild coordinates. For later comparison, let us first transform the Lemaitre 
solution from Lemaitre coordinates (r, ar) to the ordinary Painleve - Gullstrand (PG) coordinates (r, R) 
corresponding to observers moving on radial timelike geodesies labelled by Schwarzschild radial coordinate 
R. This yields a non - diagonal metric whose line element is of the form |30l [3T] 



where we used r for the Lemaitre and PG time coordinate. We present in appendix IB. 2 1 how to transform 
from Schwarzschild to PG coordinates. For a pedagogical introduction to PG coordinates and more details 
see for instance |32j . Often one uses Kruskal coordinates to obtain a maximally extended Schwarzschild 
spacetime. The PG coordinates do not cover the total Kruskal manifold. In the form given in equation 
()5.15p they cover the black hole horizon and the black hole singularity at i? = hence the physically inter- 
esting part of the manifold for our discussion. If we choose the sign of R = i\/E + F/R to be negative then 
we obtain another set of coordinates with a line element similar to that in (j5.15p but with a minus sign in 
front of the off diagonal term. These PG - coordinates will then describe the so called "white hole region" 
of the Kruskal manifold. One of the properties of Kruskal coordinates is that R is given only implicitly 





(5.15) 
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therefore for our purpose of discussing the Newtonian hmit GPG coordinates are more appropriate. One of 
the interesting properties of the ordinary PG metric shown in (jS.lSp is that the spatial hypersurfaces (that 
is dr = 0) are flat since then ds"^ = dB? + R^dQ^ and all information about the curvature is encoded in 
the shift vector given by iV = (A^^ := Rg/R, 0, 0). Furthermore, asymptotically {R oo) the PG metric 
coincides with the Minkowski metric. 



Now we transform the LTB coordinates {T,ar,0,(j)) into GPG coordinates {T,R,6,(j)). Following [29] this 
yields the metric: 



dR+ ^JE{T,R) + ^^dT^ 

ds' = -dT' + ^ — — ^ + (5^^g^ 

i + -C/(t, H) 

subject to E{t,R) > —1. The functions E and F that are functions of ar only in LTB coordinates now 
become functions depending on r and R in GPG coordinates. Therefore the time dependence of E and F 
in these coordinates is restricted by the partial differential equations 



dEir,R) 1^, , F{r,R)dE{T,R) _ 
—^-^E{T,R) + ^ - 



™-^i.(.^).^™ ^ (5.17) 

More details about the transformation from LTB to GPG coordinates can be found in appendix IB.4I We 
will not comment on the full class of solutions of this system of PDE's but just remark that e.g. (local) 
analytic solutions can be found by providing analytic initial data F{0,R), E{0,R) and determining the 
coefficients of the Taylor expansion (in terms of r) of E{t, R), F{t, R) by taking higher derivatives of (j5.17p 
at T = (Kovalevskaja method). For instance, a trivial solution is given by taking E =const. and F =const. 
In the special case E{t, R) = and F{t, R) = Rg = const the GPG metric coincides with the PG metric 
justifying the name generalised Painleve - Gullstrand metric. In contrast to the PG solution in general the 
spatial hypersurfaces of the GPG metric are no longer flat due to the 1/(1 + E) factor in front of dR^. So 
one could think that in the limit of vanishing dust density e(t, i?) — > a transformation to Schwarzschild 
coordinates is only possible in the marginally bound case E = 0. As discussed in |29j this is not the case 
and we summarise their discussion in the following: 

Let us consider the case of vanishing dust density e{t,R) = 0. Since e{T,R) ~ F' it follows F = Rg. The 
equation for F in (15.170 is thus trivially satisfied. For each function E satisfying the equation in (15.170 there 
exists a transformation from GPG coordinates (r, R, 9, (p) to Schwarzschild coordinates (T, R, 9, (j)) given by 



9ii - 1-% ^^-^^^ 

Thus one obtains a family of coordinate transformations parametrised by the functions E. Indeed, the 
integrability condition for the system (jS.lSp leads back to the first condition in (j5.17p (notice that one has 
to write E(T^ R) = E{t{T, R),R) in order to derive it). 

Therefore by first transforming the LTB to GPG coordinates it is possible to transform to Schwarzschild 
coordinates in the limit of vanishing dust density for all (allowed) functions E. Thus we can transform to 
Schwarschild coordinates not only in the marginally bound case E = but also in the elliptic {E < 0) and 
hyperbolic {E > 0) cases provided F =const. Consequently, the Newtonian limit is correctly implemented 
for all values of E and F = const. 



14 



6 Semistatic properties 



Now we want to analyse whether the physical spacetime (M, G) under consideration is semistatic in a certain 
sense or not. We already know that the physical metric G^i, takes the form of a Lemaitre metric in the limit 
e(a"r) — > which in turn is just a Schwarzschild spacetime written in comoving coordinates (see appendix 
IB.ip . So one would expect to recover its static properties in some limit. 
Let us have a look at the Killing equation 







(6.1) 



for the observable metric G^^ and examine whether there exists a timelike Killing vector field ^. For the 
Lemaitre metric such a Killing vector field is given by = [1, 1, 0, 0]"^, at least outside the event horizon, so 
we would expect to find something similar here. We will start with the ansatz 



For the physical metric described above (j6.ip reduces to 



CR + C''R' 










From ()6.5p we get 



^fr- 



aud using (|6.3p we obtain 



6.61) can be solved for 



R RR' 



in' 



and using this in (j6.7p we see that 



R' R + RR' 



it JXtt It It 



Putting all this together we obtain a condition on the solution 

R R — RR = , 
and using the explicit form of the solution ()5.12p we see that this holds only if 

F'{ar) = . 



(6.2) 



(6.3) 
(6.4) 
(6.5) 
(6.6) 



(6.7) 



(6.8) 



(6.9) 



(6.10) 



(6.11) 



(6.12) 



15 



This means that there exists a timehke Kilhng vector field ^ only if e{ar) = 0, that is if the dust matter 
fields are vanishing. 

Nevertheless, as one can easily check, the metric is invariant under the action of the vector field 

(6.13) 

up to terms of at least 0{e{ar))- So with 

(£^-;G)^, = 0{e{ar)) (6.14) 

there is a precise sense in which the physical metric G^u can be called "semistatic" . For vanishing e{ar) this 
would be an exact symmetry and the metric would be isomorphic to the standard Schwarzschild metric. 
e((Tr) = would be inconsistent because we are real dust fields as dynamical clocks rather than ideal ones 
and vanishing e{ar) would mean that the clocks do not carry any energy at alH. Nevertheless, e{ar) is a 
free function in the theory and we can choose it arbitrarily smalo. So by tuning e(cTr) in the right way we 
can get as close to a static spacetime as we wish. 



1,-^,0,0 



7 Discussion of Singularities 

In this section we will discuss the possible singularities which might appear for the metric given in (j5.9p and 
displayed again below: 

n/2 

ds'^ = -dT'^ + - ——dar'^ + R^dn^ . (7.1) 

1 + E[ar) 

As is well known from earlier studies concerning LTB spacetimes (see for example \33\ I27j ) there appear 
two different kinds of singularities: Recalling that the dust density was given by pdust = there are 

potentially two cases when it can diverge. The first one, if R{t, ar) = 0, is called a collapse singularity in 
the literature, and the second, if R'{T,ar) = 0, is known as a shell crossing singularity. It was pointed out 
by [33] for the first time that in the latter case an appropriate decay behaviour for F' when 0"^ — > oo can 
avoid the divergence of pdust in the shell crossing case. From the line element above we can read off that 
in the case of a collapse singularity the metric components Q00 and Q^cp vanish while for a shell crossing 
singularity only Qa,.a^ is zero. 

7.1 Shell crossing singularities 

For standard LTB spacetimes with positive dust densities one can always choose the arbitrary functions 
E{ar), F[ar) in such a way that shell crossing singularities will not appear ^3]. This holds not only for the 
marginally bound case {E{ar) = 0) but also for the hyperbolic {E{ar) > 0) and elliptic (—1 < E{ar) < 0) 
case. 

In our framework the dust energy density enters with negative sign. Therefore F'[a) has negative sign. 
Furthermore we have the requirement F/ R + > 0. Consequently we arrive at slightly different conditions 
on E{ar), F{ar) in order to avoid shell crossing singularities. 

The analysis for E{ar) 7^ involves the parametric solutions for R{T,ar) shown in (jS.lOp and (jS.lip . 
A singularity discussion becomes more involved in these two cases and we will restrict ourselves to the 
marginally bound case {E[ar) = 0) here, similar results can also be obtained for the hyperbolic and elliptic 

*The dust part of the Hamiltonian constraint would vanish. 

^There are certain restrictions on e{ar) which have to be fulfilled in order to avoid singularities, see section [7] for further 
comments on this issue. 
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case. 

If = then Qa^cyr is given by 



,2 1 r . i2 



F'(ct^)((t^-t) + 2F((t,.) . (7.2) 



Considering the constraints F > 0, < 0, cr-r — r > 0, in order to avoid sheh crossing, i.e. R' ^ for all 
r, ar we obviously must have 

2F 2F 

0<-F' < or -F' > > (7.3) 

CTf. — r ar — T 

for all ar, t. Since F only depends on ar, for any given ar we can choose ar — t arbitrarily small so that for 
any choice of F the second possibility can be violated. Hence the only dynamically stable condition is 

2F 

0<-F' < (7.4) 

ar — T 

Since 0",. > o",. — r > 0, this condition is certainly implied by 

9 rp 

0<-F' < — ^ < [\n{Far^)]' (7.5) 
ar 

Since the logarithm is an increasing function, we obtain that Far'^ should be an increasing function while F 
should be a decreasing function. This leaves us with a large class of possible F, for instance F{ar) = ^^"^ 
which also fulfills F < Rg. Inserting this form of F into R and R' yields the following density for the dust 

Pdust = -o /i , If — : — vTT \ 

3 (1 + i (dr + r) j [ar - T) 

The dust density diverges at t = ar as expected because this is exactly the singularity at = that can 
not be avoided, see also discussion in the next section. If we take the limit ar ^ oo for fixed values of r 
then /Odust — > 0. Hence, pdust shows a physically reasonable behaviour. 



7.2 Collapse singularities 

The second kind of singularities are those where R = and therefore besides the dust density /Odust also the 
expression for R^ in (|5.9p becomes singular. These singularities occur for all three cases - marginally bound, 
hyperbolic and elliptic - when r — > cj^. Having collapsing dust shells in mind this is exactly the moment 
when the dust shell labelled by ar reaches the singularity. Note that close to i? = the expression for R^ 
is dominated by the E = regime and thus we will again restrict our discussion to the marginally bound 
case. Looking at the explicit form of R{t, ar) written down in ()5.12p we can analyse how the different sign 
of the dust alters the behaviour during the collapse. In contrast to the shell crossing singularity there is no 
way to get around this collapse singularity for the system gravity plus pressureless dust. As long as we do 
not introduce additional matter fields and take their non-gravitational interactions into account! ^'^I we end 
up with a singularity after r = cj^, just as in the gauge variant framework using standard LTB - models. 



Finally we want to summarise all the physical selection criteria for e{ar) in the marginally bound case 

^"in phenomenological matter models this is usually done by introducing matter with non vanishing pressure. This pressure 
is mandatory in order to describe stable configurations such as stars. In the gauge invariant formalism discussed in this article 
one would essentially follow the same steps, but one has to take into account the additional component pdust when setting up 
the equation of state. For non-collapsing configurations pdust could be tuned arbitrarily small. 
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we have encountered so far: First, due to the special role of the dust fields in our formalism, we demand 
F' < or equivalently e > 0. Second we need to ensure that < F{ar) < Rs, otherwise ()5.12p is not a 

oo 

solution anymore; in terms of e this means that / d(Tre{(Tr) < 2Rs. Third, we need to choose e{ar) such that 



|-F'(iTr.)| = —F'(ar) < ^'^^'^''^ is fulfilled in order to not run into any shell crossing singularities. Given these 
conditions the physical metric G^i, is well defined until r — > cr^ when a free falling observer approaches the 
central singularity. 

8 Conclusions 

The task of finding Dirac observables for General Relativity involves the construction of quantities that 
commute with all constraints of General Relativity. Dirac observables are thus the associated to gauge 
invariant objects for General Relativity where the gauge group is closely related to Diff(A^), the group of 
diffeomorphism of the underlying manifold M. The reason why it is more complicated to construct gauge 
invariant quantities in the framework of General Relativity is that the mathematical structure of Diff (A^) is 
richer than the structures of the gauge groups used for instance in the standard model of particle physics. In 
the context of the Relational formalism one can at least formally construct Dirac observables , with respect 
to chosen clocks, one for each occurring constraint. Choosing these clocks means choosing an observer that is 
dynamically coupled to the system. Considering the case of vacuum gravity these clocks must necessarily be 
four components of the four metric g^^. If one applies the techniques developed in the relational framework 
and tries to compute (Dirac) observables for pure GR then one realises that the physical time evolution 
of those observables does not resemble the Einstein equations. This is because their dynamics is described 
by an observer who sits in a laboratory whose motion through spacetime is defined through the dynamics 
of these four metric components. Of course, theoretically there always exist a coordinate transformation 
from this observer to for instance a free falling observer, however, practically such a transformation will 
be hard to find. Fortunately, for certain types of matter coupled to gravity, the constraints, consisting of 
the gravitational and the matter contribution, can be rewritten in (partially) deparametrised form. Hence, 
these enlarged systems fall into the class of deparametrisable theories for which the question of observables 
can be addressed technically easier then in the general case. This idea was used in [19\ [20] where instead 
of considering pure gravity (i.e. the Einstein-Hilbert action) one considers gravity plus pressureless dust- 
matter fields (i.e. Einstein-Hilbert plus dust action). The dust fields become the clocks of the system and 
correspond to a dynamically coupled free falling observer. One obtains physical equations of motion for the 
observable 3-metric Qj^ and observable momenta P^^ in the dynamical reference frame defined through 
the dust fields. Of course, when introducing a dynamically coupled observer, we have to ensure that the 
occurring fingerprints of the observer are still in agreement with experimental data. 

That this is the case for the cosmological sector was already shown in |19[ [20] and in this article we 
demonstrated that this framework also describes gravitational physics in the spherically symmetric sector 
to arbitrary precision. We showed that within this sector the dynamical evolution of Qij , P*-^ , generated by 
a physical Hamiltonian Hphys, is in one — to-one correspondence with the class of Lemaitre-Tolman-Bondi 
(LTB) solutions for standard GR. Interpreted in our language, the choice of lapse and shift considered for 
the LTB solutions correponds to a gauge fixing of the (spherically symmetric) spacetime diffeomorphism 
invariance whose gauge invariant extension is precisely induced by our choice of clocks. In addition, while in 
the usual LTB framework one starts with a spherically symmetric, pressure free, perfect fiuid Ansatz for the 
energy momentum tensor whose dynamics is then derived from the Bianchi identity, here we start with a 
fully covariant matter Lagrangian whose equations of motion then reduce to the usual ones. In other words, 
our framework provides a Lagrangian underpinning of the usual LTB framework. That this works so well is 
non trivial as we discuss in appendix [Cl 

The LTB class of solutions has been carefully investigated in the literature and its dynamics in the 
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comoving frame is well understood. One might object that that the LTB class does not contain static 
solutions, at least for non-vanishing dust energy which is mandatory if one wants to use the dust fields as 
clocks. But this had to be expected when taking into account the influence of a realistic observer (i.e. an 
observer which is dynamically coupled to the system as opposed to a mere test observer). In a sense one has 
to abandon the idealisation of a static (vacuunr^) spacetime when describing physics in terms of observable 
quantities. 

Nevertheless, we showed that there exists a well defined notion of a semistatic spacetime and one can 
get as close to the standard Schwarzschild solution as desired by appropriately choosing certain constants of 
motion (which are related to the dust energy density). Additionally, we discussed that by transforming the 
LTB system to generalised Painleve - Gullstrand coordinates and considering the limit of vanishing dust 
density there exists a well defined coordinate transformation to Schwarzschild coordinates for the elliptic, 
hyperbolic and marginally bound case. Consequently, in the limit where the dust density can be neglected 
the Newtonian limit is correctly implemented. Finally, although we consider phantom dust rather than usual 
dust, there exists a range of solutions for which shell crossing singularities are avoided while the collapse 
singularity is unavoidable. 

To conclude, the framework presented in [191 [20] seems to be compatible with observations, at least in 
the cosmological and spherically symmetric sectors which are the most relevant analytically solvable ones 
when it comes to phenomenological applications. This framework might also be useful when it comes to 
quantising General Relativity. The fact that the constraints have already been solved at the classical level 
opens the door for a reduced phase space quantisation as opposed to the Dirac programme which is usually 
employed in Loop Quantum Gravity [22\ I23j. First steps into this direction have already been performed in 
|21j and a detailed analysis of this framework in the context of spherical symmetry is the subject of future 
work pi] . 
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A Boundary conditions 

We are dealing with asymptotically flat spacetimes and hence must impose suitable boudary conditions. 
For the full theory, these boundary conditions for the gauge invariant observables as obtained via Brown 
- Kuchaf dust reduction were discussed extensively in [19]. Here we need their reduction to spherical 
symmetry. Consider, as in the main text, a spherically symmetric coordinate system as (fXr, ag^a^). Tensor 
indices on the dust manifold S are denoted by i, j. A:, • • • = 1, 2, 3. The reduced observables are Qij{cr) and 
their canonically conjugate momenta P^^{a), both of which are observables in the sense described above. As 
usual, spherical symmetry constrains these flelds to the following non vanishing components and coordinate 

^^These considerations do not necessarily hold for more realistic models with additional matter fields, i.e. equilibrium states 
for complex systems such as stars still exist, just the point of equilibrium will be slightly shifted due to the influence of the dust. 
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dependence respectively 



Qij ('^r) 



diag 



(fJr ) , R^{cFr), R^icTr) sin^ 0-61 



sin (Tgdiag 



PA{ar) PRidr) PnicTr) . „2 

sm ag , 



2A{ar) ' 4:R{ar) ' 4R{ar) 



(A.l) 
(A.2) 



Following [35], in particular the parity conditions derived there, the following decay behaviour is sufficient 
to guarantee a well defined symplectic sytructure 



A 

R 

Pa 
Pr 



(A.3) 



In order to derive well defined equations of motion one should also make sure that the Hamiltonian is 
finite and functionally differ entiable. In contrast to the gauge variant framework, it is not the canonical 
Hamiltonian (with independent lapse and shift fields) but rather the physical Hamiltonian (with prescribed, 
field dependent lapse and shift) that one has to consider. While their variations are algebraically almost 
identical, the field dependence of lapse and shift prevents one from adding the usual ADM counterterms 
under the usual decay behaviour. 

In [19] it was already emphasised that, due to the dynamical nature of lapse and shift in this framework 
(see (j2.24p ). the (geometry parts of the) diffeomorphism constraints Cj must fall off strictly faster than the 
(geometry parts of the) Hamiltonian constraint C in order to accommodate asymptotically fiat spacetimes: 



lim ^ 







(A.4) 



So one has to be careful in choosing elementary fields Qij,P^^ such that these conditions hold. As we have 
seen in the last section the standard conditions are generally not strong enough for this purpose, they only 
guarantee that C — > 0{ar~^^^''^) and Ccr^ 0{ar~^^~^''^)- In this work we did satisfy this stronger fall 
off behaviour by demanding that the shift (or equivalently the momentum density) = —Q^^Cj/H = 
vanishes. It is then sufficient to add to the physical Hamiltonian the ADM mass term 



Eadm 



— Sir lim 

K r-+oo 



(Jr 



R 

h?(Tr + 2RR' 



(A.5) 



In the general case one has to choose dynamical fields A, i?, P\, Ppt such that, in addition to (|A.3p . also 
()A.4p holds. This means that the physical Hamiltonian decays as — > 0(0"^"'-^'''^^) and asymptotically 
lapse and shift behave as 



N 



H 



H 



(A.6) 
(A.7) 



B Spherically Symmetric Coordinate Systems 

In the main text we have worked with various presentations of spherically symmetric metrics in various 
coordinate systems. For the benefit of the reader we recall here how these coordinate systems are related with 
each other. Our notation is as follows: We call (r, cr^) the Lemaitre time and radial coordinate which coincide 
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with our dust time and radial coordinate. Schwarzschild coordinates are denoted by (T, R). The (generahsed) 
Painleve ~ GuUstrand hybrid coordinates are {t,R). The Lemaitre, Schwarzschild and strict Painleve - 
Gullstrand solutions are nothing else than coordinate transformations of the static vacuum Schwarzschild 
solution into comoving coordinates (on restricted patches of the fully extended Kruskal spacetime). The 
LTB family are not vacuum solutions and are expressed most easily in Lemaitre coordinates. Transforming 
the non trivial LTB solutions (i.e. non vanishing dust energy density) into Schwarzschild coordinates is not 
possible without picking up a non vanishing shift since these spacetimes are neither stationary nor static. 
However, there is a notion of semi staticity as elaborated on in the main text. 



B.l Lemaitre solution 



Starting from the usual Schwarzschild solution for vacuum spacetimes one can perform a coordinate trans- 
formation into comoving coordinates and arrive at what is known as the Lemaitre solution |28j . 

Let us start with the Schwarzschild solution in a spherical coordinate chart (T, R,6,(j)), where T and R 
approach the usual Minkowskian temporal and radial coordinates of an observer located at spatial infinity. 
In this coordinate system the line element can be written as 



l-^]dT' + 



1 



1 _ 
^ R 



dR^ + R^dn^ 



(B.l) 



where Rg = 2MG/c^ is the Schwarzschild radius, M the central mass and dQ'^ := d9^ + sin^ 9 dcj)"^ the area 
element on the unit sphere. That means the metric components are given by 
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diag 
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Then we perform the following coordinate transformation 



dT 
dar 



dT + d% 



R\l-^ 



dT + 




Rs^l 



Rs, 
R 



R^, R^sint 



dR 
dR 



or in integrated form 
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1- J# 
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2 IR 
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(B.2) 

(B.3) 
(B.4) 

(B.5) 
(B.6) 



Obviously this coordinate transformation is only valid for R ^ Rg where the Schwarzschild coordinates are 
not well defined. The inverse coordinate transformation is given by 

2/3 



R{t, ar 



3/2y/Rs{ar-T) 



T{T,ar) = T-2y/RsR{T,ar)+log(l + 

or in differential form 



IR{t, ar 



Ra 



log 



' R{t, ar 



dT 



dR 



1 



El. 

R 

R. 



-dr + 



Rs 
R 



■dar 




dr + 




dar 



(B.7) 
(B.8) 

(B.9) 
(B.IO) 
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This leads to the following line element for a vacuum spacetime in the coordinate chart (r, 



R{t, ar 



-dar^ + R (r, ar)d^} 



(B.ll) 



This solution is known as Lemaitre metric and as one can easily read of the metric components in this 
coordinate chart 



9, 



Lem 



diag 



1, 



R{t, ar 



R^{T,ar), R^{T,ar)sm^e 



(B.12) 



This metric describes the local coordinate system of an observer who is freely falling under the the influence 
of the central mass. It it not obvious (as for the Schwarzschild solution) right from the beginning that this 
solution describes a static spacetime, but one can easily calculate that for R > Rg there exists a timelike, 
hypersurface-orthogonal Killing vector ^ oc [1, 1,0,0]"^ . 



B.2 Painleve Gullstrand - Solution 



In order to transform the Schwarzschild solution shown in equation (jB.ip and ()B.2p respectively to Painleve 

- Gullstrand coordinates, we introduce an observer that moves along ingoing radial, timelike geodesies of 
the Schwarschild spacetime. Thus the observer's time is identical with the Lemaitre time r introduced in the 
last section in equation ()B.3P in differential form and in equation (jB.SP in integrated form. In contrast to the 
Lemaitre metric the Painleve - Gullstrand solution [30^1311 has as the radial component the Schwarzschild R. 
Hence, we want to perform a transformation from the Schwarzschild coordinates (T, R, 9, (p) to the Painleve 

- Gullstrand coordinates {t,R,0,(/)). This coordinate transformation has the following differential form 



dr := dT + \ 



1 



(B.13) 



R\1-^J 

The vacuum spacetime in the coordinate chart (r, R, 6, (j)) is then given by the following line element [30\ [3T] 

\ 2 



ds' 



-dr^ + { dR+ \l^dT ] + R^dn"^ 



(B.14) 



The components of the Painleve - Gullstrand metric are given by 



I 





V 



R ^ 










i?2 

R^sin^ej 



(B.15) 



This metric is no longer diagonal but still has a simple form. In particular the spatial hypersurfaces 
associated to this spacetime are flat because for dt = we get ds^ = dR^ + R'^dQ^. All information about 
the curvature is encoded in the shift vector N = {N^ = ^ Rs/ R, 0, 0). As in the case of the Lemaitre 
solution the coordinate transformation from Schwarzschild to Painleve - Gullstrand coordinates can only 
be performed when R ^ Rg. This coordinate transformation corresponds to a negative R in Schwarzschild 

coordinates (in general R = ity'^^), R decreases in time since we are considering observers moving along 

ingoing geodesies. Choosing the opposite sign for R yields the transformation 



(B.16) 
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For this form of dt we end up with the following line element 



-dr^ +\dR 




dt + R^dn^ 



(B.17) 



This line element corresponds to the so called "white hole region" of the extended Schwarzschild spacetime. 
B.3 Lemaitre— Tolman— Bondi solutions 

The Lemaitre-Tolman-Bondi solution (LTB) is a family of exact solutions to Einstein's field equations (see 
e.g. [28l [36l [25] ) that describe dynamics of a spherically symmetric spacetime filled with inhomogeneous, 
pressureless dust with energy momentum tensor T^y = poUfjUu where po is the dust's energy density and 

= [1,0,0,0]"^ its velocity vector field in comoving coordinated^. 
A general solution to this problem (in a spherical coordinate chart (r, (Jr, 0, (/>)) is given by 



ds' 

Pd(t, CJr) 
R{T,ar) 



i + £j{ar ) 

2F'{ar) 



R^{T,ar)R'iT,ar) 



±jE{ar) + 



F{ar) 

R{t, (Jr] 



(B.18) 



where slash and dot denote derivatives with respect to o"r and r respectively. 

So one can characterise a particular model by choosing particular fJr-dependent functions E{(Tr) and F{ar)- 
Then the physical radius R{T,ar) is fixed up to an arbitrary function f3{ar) which characterises different 
initial conditions for i?((7,., 0). The quantity E{ar) determines the evolution of R as well as the local geometry 
and F{ar) is related to the mass inside a shell with radial label £7^, see [25] or for a more detailed explanation 
of the physical relevance of this model. 

For different values of E{ar) one obtains different solutions, which using the notation of [26] can be given 
in parametric form as follows: 



E{ar) > 0: 



R{t, ar) 
(sinh rj — 7]) 



(coshr/ — 1) 
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2E^I\armar)-T) 
F{ar) 



(B.19) 



E{ar) = 0: 



■jF{ar){P{ar 



2/3 



(B.20) 



E{ar) < 0: 



R{t, ar) 
(r/ — sin 77) 



[1 — cos r]) 



2{-E{<Jr)) 

2{-E{ar)f/^{fi{ar)-r) 
F{ar) 



^See appendix [C] for a covariant derivation of the equations of motion for this model. 



(B.21) 
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The LTB-family of solutions is widely used to describe astrophysical situations, for example it can be used 
to model the gravitational collapse of a (only gravitationally interacting) matter cloud. Furthermore it can 
be applied to cosmological models which go beyond the standard assumptions of homogeneity in ordinary 
FRW-evolution. 



B.4 Generalised Painleve — Gullstrand Solutions 



Following [29] we present in this section how the LTB metric can be transformed into a metric expressed in 
terms of generalised Painleve - Gullstrand coordinates (GPG). As for the ordinary Painleve - Gullstrand 
coordinates we want the LTB ~ time and the GPG - time to coincide and for the radial component we 
take the function i?(r, (7^) occurring in front of dQ^ in the LTB line element in equation ()B.18|) . Hence, 
the transformation from LTB - coordinates {T,ar,0,(j)) to GPG - coordinates denoted by {T,R,6,(j)) is in 
differential form given by 

dR = —dr + —dar = Rdr + R'dar (B.22) 

OT oar 

Let us consider the following general ansatz for the line element 

ds'^ = -Xdr'^ + YdR^ + ZdrdR + R'^dn'^ (B.23) 

where X, Y, Z are functions of r and R. Considering the explicit form of dR in equation (jB.22p and comparing 
with the LTB line element in equation (jB.lSP we obtain the following conditions for the functions X, Y and 
Z 

X - YR? - ZR = I (B.24) 
YR!- = ^ (B.25) 
ZR! + 2YRR! = (B.26) 
This system of equations has the following solutions for X, Y and Z 

where substitution for R via (D.l) is being understood. The function E{ar) occurring in the LTB solution 
becomes a function of E{t, R) when expressing in terms of r, R. The equation for R? in LTB coordinates 
given in equation (|B.18p has the following expression in the GPG coordinates 



R' = E{t,R)+ ' ^ R = ±^E{t,R)+ ' . (B.28) 

Since we want to cover the black hole region of the extended Schwarschild spacetime, we choose similar to 
the case of the ordinary GP coordinates a negative square root for R and obtain the following line element 

'dR + ^E{T,R) + ^^dT\ 

ds^ = -dr^ + ^ ^ ^ + R^dn^ (B.29) 

1 + E{t, R) 

This line element is well defined for all values of r, R for which E{t, R) > —1 and E{t, R) + F(t, R) / R > 0. 
The opposite sign choice for R leads to a minus sign in front of the dt - term in the bracket of the second 
term and describes, as before, the generalised "white hole region" of the spacetime. 
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Now when solving the ADM equations that lead to the GPG metric one obtains an evolution equation for 
the shift vector N{t,R) = {N^{t, R),0,0) given by 



Cn \^{N^{t, R))^ - j = with n" = (1, -N^, 0, 0) (B.30) 

being the unit vector normal to the spatial hypersurfaces of that spacetime. This is a second order equa- 
tion for because as shown in |29j the term F{t, R)/ R is related to the shift vector component by 
2RCr-i{N^){T,R) = F{t,R)/R. See [29] for more details. From equation ([B:291) we can easily read off the 
explicit form of the shift vector 



N^{t,R) = ^E{t,R) + (B.31) 
Consequently, the equation for the lapse carries over to an equation for E{t, R) given by 



Furthermore, by applying the Lie derivative £^ onto {N^Y ^'^d using its explicit expression in terms of 
F and R shown in equation ()B.3ip we obtain 

Cn{F) = iV^ (2i?£,,(iV«)(r, R) - ^^%^) = (B.33) 

which yields a partial differential equation for -F(r, R) 



Consequently, for the GPG metric only those functions E and F are allowed that satisfy the partial differ- 
ential equations in ()B.32|) and ()B.34p . 



C Covariant analysis of spherically symmetric gravity plus pressureless 
dust 

For the sake of completeness, we analyse the dynamics of pressureless dust-matter coupled to gravity in 
a spherically symmetric setting using the usual covariant framework of GR. The metric g^^ is a priori not 
a gauge invariant object and the "dynamics" generated by Einstein's equations has to be interpreted as 
gauge transformations in the strict sense. In order to make the framework physically meaningful, we must 
fix the spacetime diffeomorphism freedom. We therefore fix a comoving coordinate system with respect to 
which the lapse is unity and the shift vanishes. In the Hamiltonian language, this completely fixes the gauge 
freedom generated by the radial diffeomorphism and Hamiltonian constraints respectively. The comoving 
coordinates will be denoted by (t, r) in order to emphasise that they are measuring proper time t along 
ideal test observer geodesies labelled by r. In the usual LTB framework, and in contrast to the Brown 
- Kuchaf framework, there is no "observer Lagrangian" that actually models these observers and their 
graviational interaction and whose proper time and geodesic label we have denoted (r, ar) throughout the 
text. Rather, one constructs a spherically symmetric, pressure free, perfect fluid energy momentum tensor 
whose Lagrangian origin remains obscure and whose dynamics is simply induced by the Bianchi identity of 
the Einstein equations. 

In what follows we will see that one ends up with an exact mathematical match between the two frame- 
works although their conceptual starting points are quite different, upon identifying (r, ar) '■= (t, r). At first 
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sight this may seem mathematically not too surprising because the energy momentum tensor of the Brown 
~ Kuchaf Lagrangian has a perfect fluid form whose pressure is constrainecj^ to vanish. However, since the 
Brown - Kuchaf Lagrangian involves altogether eight dust fields to begin with and displays a complicated 
gauge symmetry involving first and second class constraints, it is after all not straightforward to see that 
one obtains a perfect match. In particular, the velocity field of the Brown - Kuchaf Lagrangian is a com- 
plicated aggregate composed out of dust fields and a priori cannot be prescribed to take a distinguished form. 



We will assume spherical symmetry, therefore we can make the following ansatz for the line-element ds^: 

ds^ = -dt^ + K^{t,r)dr'^ + R^{t,r)d6'^ + R^{t,r)sm^ed4? . (C.l) 
This means the metric components are given by: 



ff^^ = diag -1, A\t,r), R'^{t,r), R\t,r)sm'^ 



Now we compute the Christoffel-symbols F^j^ 

abbreviation Qf^i^^^^ := -^g^Mu- 

The Christoffel-symbols are given by: 
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of gfj.u where we used the 
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(C.3) 
(C.4) 

(C.5) 
(C.6) 
(C.7) 
(C.8) 
(C.9) 



This leads to the following non-vanishing components of the Ricci-curvature tensor 



^That is, the vanishing of the pressure is an equation of motion and not put in by hand. 
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The Ricci-scalar R :- 
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Now we want to couple matter to the system, more precisely we will use inhomogeneous pressureless dust- 
matter. Assuming spherical symmetry the pressure free dust stress-energy tensor is given by 



T, 



pv 



PYi{t,r)UpUu 



(C.12) 



where pi){t,r) can be interpreted as the dust energy density and is the dust velocity vector field which 
in comoving coordinates takes the form = 5^ . 

Having collected all ingredients we can now write down Einstein's equations i?^^ — 1/2^^^ Rg^^i, = k/2T^i, 
with K = IQttG/c^ for the system under consideration: 
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and (/>(/)-comp. 
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Imposing the Bianchi identity (energy momentum conservation) we find 

y^T^u = -5VD^[ln(pDAii2)] = 



(C.13) 



(C.14) 



which obviously constrains po to have the form —e{r)/[A{t,r)R?'{t,r)\ for some free function e of r only. 

One can easily see that this system of partial differential equations coincides formally with the one 
obtained in section HI equations (j4.9p - (|4.12p . So we can map the problem of finding spherically symmetric 
solutions in the framework of gauge invariant observables to the problem of finding spherically symmetric 
solutions for gravity coupled to pressureless dust-matter in the usual framework of Einstein's General 
Relativity in the comoving gauge iV = 1, A^*" = 0. 

The system of solutions to this system of partial differential equations is the so called 
Lemaitre-Tolman-Bondi family (see appendix IB. 3|) . 
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